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Evolution of mixing of granular solids in a slowly rotated 2D drum is 

considered as a discrete mapping. The rotation is around the axis of the 

. upright drum which is filled partially, and the mixing occurs only at a free 

surface of a material. The most simple cases that demonstrate clearly the 

essence of such a type of mixing are studied analytically. We calculate the 

characteristic time of the mixing and the distribution of the mixed material 
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. over the drum. 
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O ' PACS numbers: 64.75, 46.10 

By the term "avalanche mixing" one calls the mixing of granular solids that proceeds 

00 ' 

only at a free surface layer. Mixing occurs in a slowly rotating upright drum (we shall 
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consider only a two dimensional drum) which is filled partially (see fig. 1). Granules mix 
in avalanches falling from the upper half of a free surface to the lower half. In spite of 
the simplicity of the process, nature of the avalanche mixing is nontrivial and experimental 
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mixing patterns look strikinga Nevertheless, nature of the avalanche mixing can be under- 
stood using purely geometrical approachEHi. Note that though mixing and segregation in 
rotated drums is studied intensivelylM! , the avalanche mixing is of special interest, because 
of clearness of the problem formulation. 

If one assumes that the difference between the angles of repose and marginal stability 
is little (the case of the finite difference is studied thoroughly by B. A. Peratt and J. A. 
Yorkei'i) it is possible to describe the mixing analyticallyllfl. (Note that small granules 
of the different fractions are distinguished only by color.) The fastest possible avalanche 
mixing takes place if fractions mix completely in avalanches, so the granular material is 
always mixed at the lower half of a free surface. The experiment^ appears to be surprisingly 
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close to this extremal regimen Below we shall consider an opposite extremal situation in 
which granules, while flowing along a free surface are not mixed at all. One may ask: is 
there any mixing in such a situation at all? We shall give a positive answer. The mixing 
exists (apart from a center part of the drum, of course), and it is rather efficient. Indeed, 
the fractions will be finally mixed! We shall describe dynamics of such a mixing. 




FIG. 1. The avalanche mixing scheme. For an infinitely small turn of the drum, the granules 
of different fractions flow from sector A to sector B, undergoing mixing. The free surface of the 
material is tilted at the angle of repose all the time. The pure fractions are shown by the black and 
white colours. Regions of mixed material are denoted by the grey colour. The degree of mixing 
is not shown. The material inside of the dashed circle will be never mixed. The result of the 
avalanche mixing depends on the way in which the granules flow from sector A to sector B. 

Then there are two possible cases. In the first one, granules with the linear coordinate 
x at the upper half of a free surface (we assume that at the free surface midpoint x — 0) 
fall to the point x' — x at the lower half, so in fact the surface wedge is not changed after 
avalanches (see Fig. 
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FIG. 2. Limit regimes of the evolution of granules while flowing from the upper half of a free 
surface to the lower half. 

a - the material mixes completely in the avalanches; 

b - after the flow, the material appears to be inverted, and upper granules become lower ones - 
the x 1 = x transformation of the wedges (x is a linear coordinate at a free surface, counted from 
the midpoint of the free surface in both directions). 

c - the material slips as a whole from the upper half of a free surface to the lower half, so lower 
granules stay to be lower ones after the slip. Using the condition of volume conservation one 
obtains easily the x' = v sin 2 9 — x 2 transformation of the wedges. Here the angle 6 characterizes 
the relative volume of an empty space in the drum. 

In the second case, the material from the upper wedge slips to the lower one as a whole 
(Fig. 0,c). The drum radius is unit, and the measures of the amount of the material are the 
angles 6 < n/2 for more then half filled drum and d = -k — 6 < h/2 for less then half filling 
(see Fig. |l|). Thus, the coordinates of the lower and the upper points of a free surface will be 
sin# or sin^. Then, accounting for volume conservation, one describes the transformation 
of the coordinates as x' = V sin 2 6 — x 2 (or x' = V sin 2 $ — x 2 ) in these situations. 
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Introducing the scaled coordinate z = x/ sin8 or z = xj sin-$ we get for the cases under 
consideration: z' = z and z' = \Jl — z 2 . 

Let us describe by the number k events of appearing of some granule at the lower half 
of a free surface. The corresponding coordinate of the granule at the lower part of a free 
surface and time will be (xk,tk). (The angle of the drum turn plays the role of the time.) 
Transformation of this pair describes completely the mixing, since all the granules move as 
a whole in the bulk of the material. We shall study this discrete mapping. 

We start with consideration of the z' = z transformation case for less then half filled 
drum. One can see from Fig. |l] that the coordinate is unchanged during the mapping and 
tfc+i = tk + 2 arctan(xfc/ cosfi). 




FIG. 3. The avalanche mixing in the event of a less then half filled drum and the x' = x 
transformation of the wedges. A free surface is shown to be horizontal, since our results do not 
depend on the angle of repose. The angle $ < ir/2 characterizes the relative volume of the granular 
material. If a granule is at the point x at the left half of a free surface, it will always stay at a 
dashed circle sector. 

Thus, one obtains the following trivial discrete mapping^: 

(zk+i, h+i) = (zk, tk + 2arctan(z fc tan$)) . (1) 

One may study the simplest case of an infinitely thin layer of the black fraction over the 
white fraction at the initial moment, that is (zo,to) = (-2,0). Then the black material will 
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evolve in the following manner: 



{zk-, tk) = (z, 2k arctan(ztant?)) . 

Thus, the black granules will be at the point z of a free surface at moments tf.(z), k 
0,1,2... (see Fig. @). 

25r 
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FIG. 4. Evolution of a thin layer of the black fraction for a less than half filled drum in the 
case of the z' = z transformation (z = /xsiwd). i? = 0.4-7T. (See Eq. (0)). Crossings of the line 
i = const with the curves correspond to the points at the lower half of a free surface in which black 
granules appears at the time t. 

Formally speaking, there will be no true mixing of the material at any time, but near the 
fixed z the distance between black granule layers are diminished with time. When, in the 
vicinity of z, this distance becomes to be of the order of the granule diameter, the material 
is in fact mixed in this area. A very simple geometrical consideration gives for this distance: 
2 zsin$(l + z 2 tan t?) arctan 2 (,ztan??) 



d 



(3) 



t { [1 + (1 + z 2 tan 2 ) arctan(^ tan t?) /tf + z 2 tan 2 ■&} 1 / 2 ' 
If now d is the relative diameter of granules (recall that the drum radius equals 1) then one 
obtains the corresponding characteristic time 
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4z 2 sin 2 $ 



1 1/2 



1 + 2 2 tant9) arctan^ztant?) — z 2 tan 2 i9 
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The expression is most simple in the following limit cases. If ztawd 1, 

d 
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If zt&ntf > 1, 



tmix 



7l 2 Z 2 
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tan sin 
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Finally, the black fraction from the point z will be distributed homogeneously over the 
dotted line shown in Fig. [| If initially the black fraction occupies the sector of the 7^0 
angle at the lower half of a free surface, then the concentration of the black fraction in the 
mixture will be 



c(z 
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2 arctan^tant?) 

so the mixing differs greatly for different z. As it can be seen from Eq. 
t m ix ~ > oo, and any mixing is practically absent (see Fig. f|). 
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when d — > tt/2, 
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FIG. 5. Same as Fig. for a more than half filled drum in the case of the z' = z transformation 



(z = sin ( 



0.4tt. (See Eq. 



The case of a more then half filled drum 9 = tt/2 — d < tt/2 can be considered in the 
same manner. In the event z' = z one writes 



tk = 2k[ir — arctan(ztan#)] 



(8) 



(see Fig. |5]) - cf Eq. (0). The material in the point z will be mixed to a "homogeneous" 
state with the black fraction consentration 



c(z) 
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2tx — 2 arctan(z tan l 



(9) 



in a time 



tmix = [tt — arctan(ztan( 



Az 2 sin 2 9 
I 2 



1 + z 2 tan 2 6) [n — arctan(z tan 8)} 2 — z 2 tan 2 6 



1/2 
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Thus, if ztanO ^> 1, 



tmix 



71 



— tan 9 sin 6 
2d 
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and in the event of ztan6 <C 1 



_ 2n 2 z . 

tmix , Sin u 
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Note, that Eqs (0) and (0) are valid only if there are many crossings of the black fraction 
trace with the lower half of a free surface, i.e. if t m ix 3> 2ti{2tt /29). Thus, as it follows from 

- 27i 2 /e. 



Eq. (|T2"D, the 6 angle should be 6 ^> Jd/z. Otherwise t mi 

We see that there is a minimum of t mix at some intermediate filling value. A full de- 
pendence of t mix on i? (or rather t m \ x ($)) is shown in Fig. ||. It somewhat resembles the 
characteristic mixing time behavior in the case of a full mixing in the wedgesS, though the 
meaning of t m i X is quite different now. 
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FIG. 6. t^] x vs i9/tt. The later one characterizes the relative volume of the granular material: 
the drum is empty, half rilled, or full when # = 0, vr/2, or tt, correspondingly. The crossover region 
near i? ~ 1 is shown by guess-work, z = 0.5, d = 0.01. 
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Now we proceed with the case of a less then half filled drum and the z' = \/T 



transformation. One see that 



(z k +i, tk+i) = (j/l - z%, t k + 2arctan(z fc tan$)^ , 



(13) 



so if initially (z ,t ) = (z,0), then 



(z 2n , t 2n ) = \z, 2narctan(zfctan-(9) + 2narctan(yl — tan$) 
(z2n+i, t2n+i) = ( V 1 — £ 2 , 2(n + 1) arctan(,2 fc tan ■#) + 2narctan(yl — ^tan??) ) , (14) 



where n = 0, 1, 2, . . .. Therefore, 



t2n = 2n arctan^tan-i?) + 2narctan(\/l — z 2 ta.n$) 



hn+i = 2narctan(ztan , i9) + {2n + 1) arctan(\/l — z 2 tani?) 



(15) 
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FIG. 7. Same as Figs || and || for a less than half filled drum and the z 1 = \J\ — z 2 transformation 
of the wedges. 9 = 0.4vr. (See Eq. pi)). 



Now there is a point z — 1/ \[2 in which the fractions will be never mixed, and mixing 
in fact occurs in the vicinity of z = or z = 1. The region of poor mixing is broaden as the 
drum filling tends to the one half level. The averaged concentration of the black fraction is 



c(z) 
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arctan(ztan'i9) + arctan(-\/l — z 2 tan$) 



(16) 



(cfEq. (0)). 



In the same manner can be considered case of the z' = \/l — z 2 transformation for a 
more then half filled drum. One can obtain easily 



t2n = 4n7r — 2narctan(ztan#) — 2narctan(vl — z 2 ta.n8) 

^2n+i = (2n + l)27r — 2narctan(ztan^) — (2n + 1) arctan(Vl — z 2 tan#) (17) 




FIG. 8. Same as Figs ^, ||, and for a more than half filled drum and the z' = \/T" 
transformation of the wedges. 6 = OAtt. (See Eq. (JO)). 



The averaged concentration of the black fraction will be 



c(z) 
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(18) 



2tt — arctan(ztan^) — arctan(vT-^tan^) 
(cf Eq. (|9|)). Note that in the event of the z' = \J\ — z 2 transformation, mixing is far poorer 
then for the z' — z transformation. 

In conclusion, using discrete mapping approach, we considered analytically two limit 
regimes of avalanche mixing, when the mixing in the wedges is absolutely absent. Neverthe- 
less, finally, the material appears to be rather effectively mixed. The experimental situationQ 
- the mixing of salt grains - is between full mixing in wedges and what we call the z' = z 
transformation of the wedges (see the experimental patterns of the mixing in the paper of 
G. Metcalfe et at). There was a question&B: why is the experiment so close to the extremal 
(the fastest!) regime of avalanche mixing, in which granules fully mix in avalanches? Why 
do the theory for this extremal regime give such a surprising coincidence with the mixing 



time values observing in the experiment? Now we can see the reason. As we have shown 
in the present communication, the dependence of the quantity characterizing the time of 
mixing on the drum filling (see Fig. |6|) in the case of the z' = z transformation appears to 
be rather similar to the behavior of mixing time in the fastest possible for avalanche mixing 
regime. Thus the experiment, which is between this two limit cases with a similar behavior, 
can be described well. 
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